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ABSTRACT 

We present mid-infrared spectra and images from the GEMINI-N (+Michelle) 
observations of fragments SW3-[B] and SW3-[C] of the ecliptic (Jupiter Family) 
comet 73P/Schwassmann-Wachmann 3 pre-perihelion. We observed fragment B 
soon after an outburst event (between 2006 April 16 - 26 UT) and detected 
crystalline silicates. The mineralogy of both fragments was dominated by amor- 
phous carbon and amorphous pyroxene. The grain size distribution (assuming 
a Hanner modified power-law) for fragment SW3-[B] has a peak grain radius of 
a p ~ 0.5 /im, and for fragment SW3-[C], a p ~ 0.3 /im; both values larger than the 
peak grain radius of the size distribution for the dust ejected from ecliptic comet 
9P/Tempel 1 during the Deep Impact event (a p = 0.2 /xm). The silicate-to-carbon 
ratio and the silicate crystalline mass fraction for the submicron to micron-size 
portion of the grain size distribution on the nucleus of fragment SW3-[B] was 
1.341 tolss and 0.335±o;m, respectively, while on the nucleus of fragment SW3- 
[C] was 0.67llg;Q7g and 0.257±q;q43, respectively. The similarity in mineralogy 
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and grain properties between the two fragments implies that 73P/Schwassmann- 
Wachmann 3 is homogeneous in composition. The slight differences in grain size 
distribution and silicate-to-carbon ratio between the two fragments likely arises 
because SW3-[B] was actively fragmenting throughout its passage while the ac- 
tivity in SW3-[C] was primarily driven by jets. The lack of diverse mineralogy 
in the fragments SW3-[B] and SW3-[C] of 73P/Schwassmann-Wachmann 3 along 
with the relatively larger peak in the coma grain size distribution suggests the 
parent body of this comet may have formed in a region of the solar nebula with 
different environmental properties than the natal sites where comet C/1995 01 
(Hale-Bopp) and 9P/Tempel 1 nuclei aggregated. 

Subject headings: comets: general: comets: individual (73P/Schwassmann-Wachmann 3) 
:dust: modeling 



INTRODUCTION 



Comet nuclei contain some of the least thermally- and least aqueously-altered materials 
extant in the solar system and are highly porous agglomerates of ice and dust grains, perhaps 
with highly stratifi e d, inhomogeneous layers of varied density, porosity, and composition 
( IHarker et al.l l2007t iBelton et al.l 120061 ) formed by layered accretion of icy planetesimals. 
The latter may or may not have aggregated from specific radial zones in the protoplanetary 
disk to form the nascent nuclei of comets, ultimately giving rise to the two comet families, 
ecliptic comets (ECs; including Jupiter Family comets and active Centaurs) and nearly 
isotropic comets (NICs; including long-period Oort cloud and Halley-type comets). 

Recent dynamical models by the 'Nice group' indicate that Jupiter ejected bodies from 
the solar system's gravitational potential and into the Oort Cloud to become NICs. Comet 
nuclei that formed in the trans-Neptune region, in situ in the Kuiper Belt (K B), and the possi- 



Morbidel 



i et al. 



20051 ) 



ble few ejected by Jupiter in to the KB, survived dynamically to become ECs 
2004 : IWhitman et al.l 120061 ). As collisions are frequent in the Kuiper Belt (iPan et al. 
Jupiter-family comets (ECs) could be chips off of larger Kupier Belt Objects (KBOs). There- 
fore, the structure of EC nuclei may be very different from NICs. On the other hand, Port 
cloud (NIC) comets (i.e., long-period) were grav itationally s cattered to the Oort cloud ([Oort 
195ll ) where they suffered no further collisions (jSternl 120031 ) . On ce perturbed into the inner 
solar system, NICs have short life ti mes compared to EC comets (A'Hearn et al.l 1985 ). ECs 
have notably less activity than NICs (lA'Hearn et al.lll985t iMeech & Svorenll2004f ). Lower ac- 
tivity and lower dust production rates result in less thermal dust emission making dust prop- 



erties of ECs relatively unknown compared to brighter NICs ( IWooden et al.ll2007 
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2004 ; lHanner. Lynch, fc Russell! Il994f ). 



The strongest tests for new dynamical models of solar system formation processes are 
the compositions of comets, and in particular, the diversity of the gas-to-dust ratio and 
dust composition within any given comet's coma coupled with variations in the latter phys- 
ical characteristics (among others) between comets within a dynamical family and between 
dyna mical famil i es. Notably, solar nebula models for thermal processing and radial trans- 
port (IBossi 12004 ; ICuzzi et al.lll993l) of grains are constrained by th e crystalline silicate mass 
fraction deduced for comets (IWooden et al.l 120071 ; IWoodenl 120051 ). as well as the (contro- 
versial) presence of aqueous altered grain components. Amorphous silicates in co mets are 
thought to be presolar ISM silicate g rains (|Woodenl 12002c lHanner &; Bradelyl 120041 ) because 
of their high abundance in the ISM ( ILi fc Draindl200ll ) and ubiquitous pre s ence as glassy 



spherules (GEMS) in cometar y interplanetary dust particles (lBradleylll994t iBradley et al. 



19991 ; lHanner fc Bradelyl I20041). On the other hand, crystalline silicates are relatively rare 
(<2.2 - 5%) in the ISM (kemper et all liooi l2005[ ILi fe Drainel boOlh yet are abundant 
in many co mets. Crystalline silicates in comets are either high-temperature (^1450 K) 
condensates (iGrossmanl Il972l). or amorphous silicates ann ealed into crystals at ~1000 K 



flHallenbeck et al.lll998l : IWooden 



<0.3 Myr-old) solar nebu la (iBockelee-Morvan et al 



200 



5; 



Wooden et al. 



20071) in the inner zones of the early 
2ooi), or annealed in nebular shocks 



at 5-10 AU ( < 1 Myr-old) (jHarker fe Deschl 120021 ) . In models of turbulent mixing, the early 
formation of crystals in the hot inner zones of the solar nebula are concurrent with rapid 
disk spr eading (0.3 Myr) suggests a u niform crystalline-to-amorphous silicate ratio between 
comets flBockelee-Morvan et al.l 120021 ) . In contrast, viscous dissipation models may yield a 
radial gradient in the crystalline silicate mass fraction that spreads to only 20 AU in 1 Myr 
(IWooden et al.ll2007l ; ICieslall2007l ). making it harder to explain crystals in comets that formed 
in the trans-Neptune zone. Strong constraints for nebula models include the high crystalline 
silicate mass fracti on (~25%-50%) deduced f rom modeling NIC C/1995 01 (Hale-B opp) 
jHarker et alil2004h . NIC C/2001 Q4 (NEAT) (barker etaflbood IWooden et all 120041 ) . the 



ejecta from the EC 9P/Tempel 1 duri ng Deep Impact (jHarker et al. 



2005, 



dant crystalline silicates (jEbel et al. 



2007m . and the EC 



78P/Gehrels 2 ( Watanabe et al.l 2005). Sta rdust samples from EC 81P/Wild 2 contain abun- 



20061 ) and even higher temperature calcium-aluminum 



inclusion (CAI)-type minerals. In contras t, other ECs have a low (^C50%) cr ystalline silicate 



mass fractions: pre- im pact 9P/Tempel 1 (ILisse et al.ll2006l ). 1 03P /Hartley 2 (ICrovisier et al. 



2000l ; ICrovisierlll999l ). and 29P/Schwassmann-Wachmann 1 (jStansberry et al.ll2004l ). Thus, 
assessing the composition of cometary nuclei is crit i cal to understan ding the origins of comets 
(jEhrenfreund et al.l 12004c lHanner fc Bradelyl I2004J ; ICrovisierll2007l ) and the formation of the 
solar system. 



Contrasts between EC and NIC family grain properties often rely on comparisons be- 
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tween ground-based remote-sensing studies and the in situ sampling of coma dust in comets 
such as lP/Halley, 81P/Wild 2 or the 9P/Tempel 1 (henceforth 9P) Deep Impact ejecta anal- 
ysis. Frequently, the 10 /im spectra of short-period ECs exhibit a br oad and weak 10 
featu re, of order 10-20% above the 8 to 13 /zm blackbody continuum (IHanner fc Zolensky 
20101 ) and little evidence for the presence of crystalline silicates - a result of space weathering, 
thermal processing, and de-volatilization arising from frequent pe rihelion passage, whereas 
in NIC comets such emission is common (jKelley fc Wooden! 120091 ). Surprisingly, the Deep 
Impact event excavated sub-surface materials containing a significant population of small 
grains, including crystalline silicate species, unlike the grain population dominant in the 
quiescent coma of the EC 9P. The observed 10 and 20 /im spectra and composite thermal 
spectral energy distribution (SED) models of the 9P dust ejecta were similar in charac- 
teristics to those d e duced for the coma grain s in comet C/199 5 01 (Hale-Bopp), a NIC. 
Harker et al.l (120071 ). ISugita et al.l ( 120051 ). and iLisse et al.l ( 120061 ) among others conjectured 
that the interior of 9P may be comprised of less processed, more primitive, or pristine grain 
material. Therefore, inference between comet families drawn from the study of coma dust is 
a challenge. The assertion that ECs are compositio nally distinct from NICs as a class is no 
longer valid based on the Deep Impact experiment (IHanner fc Zolensky 1 120 10h . 



In its 1995 apparition, the Jupiter-family comet 73 P / Schwassmann- Wac hmann 3 (hence- 
forth SW3) exhibited a sudden increase in brightness (jCrovisier et al.l 120061) near peri helion 
breaking into three bright fragments and several fainter components ( iScotii et al.lll996l ). The 
observing geometry of SW3's next periheli on passage (orb ital period T = 5.36 yr) was not 



very favorable for Earth-based observer s (ISekanina 
perihelion rj 66 fragments wer e detected (jWeaver et al. 



2005). During SW3's 2006 return to 



20081 ). many of which were remnants 



of the splitting event of 1995 (IReach et al.ll2009l ). Fragments SW3-[B] and [G] continued to 
split, and swarms of mi ni-comets were disco vered in Hubble Space Telescope images of the 
tails of these fragments (jWeaver et al.l 120081 ) . Fragmentation of this comet's nucleus, a nat- 
ural analogy to the Deep Impact mission, prov ided a rare opportunity to observe potentially 
pristine material from the progenitor interior ( IDella Russo et al.l 120071 ) . 



Here we present mid-infrared narrowband imagery and spectroscopy (§2]) and coma dust 
mineralogy derived from thermal grain models (§ |3]) of the two largest surviving fragments, 
SW3-[B] and SW3-[C] during the 2006 apparition, discussing comparison (§3} between grain 
properties of this comet (one of the most active known ECs) to the putatively more primitive 
dust populating the coma of NIC comets such as C/1995 01 (Hale-Bopp). Our conclusions 
are summarized in 35] 
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2. OBSERVATIONS 

Mid-infrared (IR) 10 /zm spectra and 10 and 20 /zm images of fragments SW3-[B] and 
SW3-[C] of comet SW3 pre-perihelion (perihelion 2006 June 08 UT) on 2006 April 29.5 UT 
(fragment SW3-[B]) and 2006 April 30.5 UT (fragment SW3-[C]) were obtained using the 
Michelle imaging spectrograph on the 8-m Fredrick C. Gillett telescope (Gemini-N) on Mauna 
Kea, Hawaii conducted under program GN-2006A-DD-1. In addition, we obtained 20 fim 
spectra of fragment S W3- [B] . The observational summary is shown in Table [TJ 



Michelle has several mid-IR broad and narrow band filters (lAfistair et al.lll997l ). For our 
observations, we used two narrow band filters: 1) the Si-5 filter with a central wavelength 
of 11.6 fim and width of 9.5%; and 2) the Qa filter with a central wavelength of 18.1 jum 



and width of 10.7%. Common mid-IR observing techniques (e.g., I Joyed Il992l ). chopping 
the secondary (±15" at a 45° angle in the NW-SE direction) and nodding the telescope 
(±15" in the NW-SE direction) to reduce the thermal noise from the background and avoid 
contamination from the coma extending to the south west for both comets, were utilized. 
The brightness of the coma in the NW-SE direction for both fragments drops to the value 
of the background at 11.2 /zm (40.12 mJy arcsec -2 ) at 6" away from the centralized core. A 
final image is produced by subtracting the two nod pairs. 

Michelle uses a grating to produce spectra with wavelength coverage of ~ 7.8 — 13.2 /im 
and ~ 17.0 — 21.0 /jm within the 10 /zm and 20 /zm ground-based atmospheric windows, 
respectively. A 0.6" wide slit was used for all of our spectral observations, providing a 
resolution of R ~ 200 and R ~ 110 at 10 /xm and 20 /im, respectively. We chopped on array 
at the maximum allowable chop throw of 15" along the slit direction, perpendicular to the 
dispersion axis of Michelle. 

Extraction of the comet and standard st ar spectra fr om the two-dimensional images used 



common mid-IR data reduction techniques (jJoycelll992l ): 1) two nod pairs are subtracted to 



produce a background subtracted image which contains both a positive and negative spec- 
trum; 2) two extractions (spatial dimension of 1.0" yielding an effective aperture diameter 
of 0.87") are made to produce each spectrum: one is the source spectrum, and the other is 
extracted 6" to the north of the source spectrum for the purpose of calculating the residual 
background; 3) all of the counts at each pixel (i.e., wavelength) inside the source aperture are 
summed; 4) the sky extraction is used to calculate a median residual background value and 
variance, which is used to calculate the photometric error for the object spectrum; and 5) the 
median sky value is subtracted from the object spectrum. Nominally, a residual sky spec- 
trum to the south of the source spectrum is also used to calculate the background value and 
variance. However, in this case, for both fragments, there is emission from the extended coma 
in this southern direction. Therefore, we have relied only on the northern sky extraction. 
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Using this technique does not affect the overall value of the sky, but simply does not enable 
better constraint of the uncertainties in the sky variances. This uncertainty is reflected in 
the larger 1-sigma error in the final spectrum. We increased the signal-to-noise of the data 
with a 3-point statistically weighted average boxcar function. This operation decreased the 
spectral resolution to R = 129 at 10 p,m. The spectral resolution is suffic ient to identify 
mineral resonances that may exist (IHarker et al.l l2002t IWooden et al.lll999l ). In 100.8 s of 
on-source integration time, a given flux point in the final spectra has a root-mean-square 
(RMS) error of 2.5 x 10~ 19 W cm" 2 /im" 1 . 



Cohen standards (jCohen et al.lll996l ) w hich provid e superior IR spectrophotometric cal- 



ibration, and ATRAN atmospheric models ( jLordlll993l ) are used to flux calibrate the spectra 
and correct for atmospheric extinction. The object spectrum is multiplied by the Cohen 
stellar template, divided by the extracted instrument spectrum of the standard, and finally 
multiplied by the ratio of the atmospheric transmission spectra calculated for the airmass of 
the comet and standard star. Photometric standards were also used to assess the seeing and 
derive estimates of the point-spread function (PSF). The average PSF FWHM was 0.4" at 
11.6 jum and 0.5" at 18.1 pm. 

Figure [1] shows the 11.6 and 18.1 /im images of fragments SW3-[B] and SW3-[C], re- 
spectively. The 3-sigma detection limit of the coma of SW3-[B] extends ~ 15" from the 
central nucleus region towards the south western direction at a position angle (PA) ~ 205° 
in both images. Notably, there is a region of enhanced flux (32% above the expected flux of 
a standard 1/p decline in coma surface brightness, where p is the projected distance from the 
nucleus) approximately 3" away from the centralized condensed region that is observed at 
both wavelengths. The coma of SW3-[C] is fairly condensed (0.6" at 11.6 pm), the 3-sigma 
detection limit of the coma extending only about ~ 8" at PA~ 219°. We note that more 
sensitive Spitzer MIPS24 observations of 73P, obtained less than a week after our Michelle 
observations (2006 May 4-5 UT), showed that the surface brightness of the coma of frag - 
ment SW3-[B] extended some 4' while that of SW3-[C] extended ~ 5' faeach et al-lbopgh . 
Isophote-contour maps derived from our SW3-[C] image do not reveal any enhanced regions 
of flux density, in contrast to the surface brightness morphology evident in the coma of 
fragment SW3-[B]. 

Figure |2] shows the location and orientation of the 0.6" wide spectral slit during the 
observation of each fragment, the location and size of the extraction areas for the spectra, 
and the extracted spectra. The slits were centered on the peak of the azimuthally-averaged 
coma surface brightness profile, essentially coincident with the nucleus (unresolved at our 
platescale), and rotated to align along the direction of the extended coma. Individual spectra 
were extracted in 1.0" segments in the spatial direction for mineralogical and grain analysis 
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The spectral extraction apertures centered on the nucleus contain both nuclear and 
coma fluxes. In order to study the coma dust properties, we need to properly assess the 
nuclear contribution to the flux. We can deconvolve the relative contribution of the nucleus 
flux to the coma by analyzing the high spatial resolution images at 11.6 pm. We take the 
PSF of the standard star HD 127665 as the shape of the nuclear contribution to the images. 
We assume a standard 1/p coma surface brightness decline. We derive the azimuthally- 
averaged cumulative photometry versus increasing p, i.e., from the peak flux pixel, for the 
image and for the sum of a scaled PSF and the 1/p-coma. By x 2 -fitting the shape of the 
azimuthially-averaged coma profile, we constrain the contribution of the PSF to the image. 
By putting a synthetic beam of 0.6" x 1.0" on this solution, the nucleus is found to contribute 
1.4% ± 0.2% of the measured flux in the centered slit at 11.6 /im to a confidence level of 
99% for fragment SW3-[B]. Our calculated scaling factor corresponds to an effective radius 
of the nucleus of 0.2 km ± 0. 05 km for fragment SW3-[B], equal to the lower limit of the 



nucleus radius determined by iHowell et al.l ( 120071 ) . The nucleus contribution is subtracted 



from the central extracted spectrum of fragment SW3-[B], so we only need to model the 
remaining coma grains. From the 11.6 pm image of fragment SW3-[C], we determine that 
the flux at the central region is produced solely from the surrounding coma, and that there 
is no contribution to the flux from the embedded nucleus (at the 99% confidence level). 
Therefore, we do not apply a correction to the spectrum for fragment SW3-[C]. 



MODEL OF THERMAL EMISSION FROM DUST GRAINS 



The lHarker et al.l ( 120021 ) thermal emission dust code was used to model the observed 
dust grain emission in the coma of fragments SW3-[B] and SW3-[C]. The model assumes 
that an optically-thin collection of discrete (singular mineralogy) dust particles reside at 
the heliocentric (r^) and geocentric (A) distance of the comet at the epoch of observations 
(Table [TJ. The 10 pm mineralogy us ed in the model is d erived from labora tory studies of 



interp lanetary dust particles fIDP s. IWooden et al.l 120001) . micrometeorites (IBradley et al. 



19991 ). the N ASA Stardust mission (jBrownlee et al.ll2006l) and other remote sensing mo deling 



efforts (e.g., lHanner. Lynch, fc Russelllll994j ; lHarker et al.l 120021 : IWooden et al.l 120041 ) . Our 



selection of basic dust grain components is consistent with the major mineral groups used by 
other modelers to generate synthetic SEDs arising from dust thermal emission in the 10 pm 
region (IMin et al.ll2005l ; ISugita et al.ll2005l ; iLisse et al.ll2006l ; IWerner et al.ll2006l ) . Amorphous 
silicates with chemical composition (stiochiometry) similar to olivine (Mg^,Fe(i_ y ))2Si04 
and pyroxene (Mg x ,Fe(i__ x .))Si03 with x = y = 0.5 (i.e., Mg/(Mg + Fe) = 0.5) repro- 
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duce the broad width of the 10 /im feature. Mg-rich orthopyroxene is detected through 
its 9.3 and 10.5 /im features (jWooden et al.lll999l ; lHarker et al.ll2002l ). Mg-rich crystalline 
olivine is uniquely identified thro ugh its distinct, relatively narrow 11.2 /im silicate feature 
(IHanner. Lynch. &: Russell! 1 1994J ) . Mg-rich crystalline sp ecies are defined as grains with a 
stoichometry of 0.9 < x<1.0 {Mg/{Mg + Fe) = 0.9) JWoodenl liool iKoike et all I2OO3I : 
Chihara et aJbooj l. 



Lastly, the presence of amorphous carbon grains in the coma is required in order to repro- 
duce the underlying featureless emission (continuum) in the 8 — 13 /im wave length region. 



Othe r species such as P AHs, hydrous silicates (clays) and iron sulfides (see I Werner et al. 



20061 ; iGehrz et al.l 120071 ). cannot be uniquely identified from our SW3 spectra. Therefore, 



these latter species are not included our model. 

The amorphous grains are assumed to be spherical (Mie Theory) and range in size 
from 0.1 — 100 /im. The indices of refraction used in Mie Theory to calcu l ate th e optical 
efficiencies of each amorp hous si l icate materials are from iDorschner et al.l (119951 ) and the 
amorphous carbon is from lEdohl (119831 ). The submicron sized crystalline olivine grains are 
modeled using a continuous distribution of ellipsoids (CDE) with one of the crystallographic 
axes elongated to a rati o of 10:1:1 to matc h the wavelength of the emission features with 
laboratory experiments (IFabian et al.ll200ll ). The use of elongated grains is supported by 
polarization measurements of the diffuse interstellar medium which suggest the presence 
of elongated grains, and con densation experime nts which show the formation of elongated 
forsterite (Mg-pure) crystals (IFabian et al.l 1200 ll ). The radii of the crystalline olivine grains 
ranges from 0.1 — 1.0 /im; large r crystals do no t well match the shape of the observed 
resonances in comet Hale-Bopp (IMin et al.l 120051 ) . The indices of refraction for the three 



crysta llographic axes of the crystalline olivine between 3 and 250 /im are from iJaeger et al. 



(119981 ). The non-oriented crystalline olivine indices ranging from 0.2 — 2 /im are from the Jena 
Laboratory Group0. The short wavelength indices are appended to each long wavelength 
axes prior to the CDE calculations. Note that the radius of the crystals is an "effective" 
radius, or the radius a crystal would be if it were spherical. 

The temperature of each grain is calculated using a t hermal equilib ration calculation. 
In the coma, the grains reach radiative equilibrium quickly (lWoodenlll979l ) as they enter into 
the coma. The temperatur e of the cryst a lline o livine grains are calculated using the "hot 
crystal model" presented in lHarker et al.l (120021 ): the temperature of the crystalline olivine 
grains are increased by 1.9 times over those predicted from the optical constants in order 
to fit the Infrared Space Observatory (ISO) SEDs from comet C/1995 01 (Hale-Bopp) at 



1 http://www. astro. uni-jena.de/Laboratory/Database/databases. html 



-9 - 



2.9 AU (jCroviser et al.l 119971 ). The temperature increase was required to match the relative 
emission peak strengths in the 10 and 20 /im regions, and may be caused by crystals being 
in contact with warm material, or the crystals themselves contain some Fe making them 
Mg-rich, but not Mg-pure. We retain this "hot crystal model" in all of our comet modeling 
to compare results between comets, although the model awaits furthe r constraining wit h 
additional long wavelength data from the comets observed by Spitzer ( iKelley et al.ll2006l ). 
An effect of using warm crystals in thermal models is the calculated dust species masses 
derived from mineral decomposition analysis are less than the mass fractions computed if 
cooler temperature grain populations were invoked. However, use of warm crystals does not 
affect the relative amount of crystalline silicates we calculate in the fragment comae of SW3. 

For the grain size distributi on, we use the Hanner modified power law (hereafter HGSD; 
Hanner. Lynch, fe Russell! 1 19941 ). In the HGSD, the relative number of each grain size is 
calculated by: n(a) = (1 — a /a) M (a /a) N , where a is the grain radius, a = 0.1 /im, the 
minimum grain radius, and M and iV are independent parameters. Since N weights more 
of the larger radii particles and can be considered as the size distribution slope for large 
radii, it is best constrained at wavelengths longer than 13 /im. The value for N is set at 
three different values in the model: 1) N = 3.4 which weights more of the larger particles in 
the size distribution; 2) N = 3.7 is the canonical va lue we determined from our modeling of 
comet C/1995 01 (Hale-Bopp) (IHarker et al.ll2002l ); and 3) N = 4.2 which weights more of 
the smaller particles in the size distribution. From the values of N and M, we also calculate 
the peak of the HGSD: a p = (M + N)/N. The value a p is not a freely varying parameter, 
but a calculated property of the size distribution. We varied M such that a p changed by 
increments of 0.1 jum. 

The amorphous dust can be modeled as porous grains (the crystals always remain solid). 
The fractal poro sity of the amorphous grains is varied to create porous or "fluffy" particles 
( Lisse et al.lll998l ). The amount of vacuum is: / = 1 — (a/ a o y D ~^ with D = 3.0 correspond- 
ing to solid grains and D = 2.5 corresponding t o the maximum fracta l porosity allowed while 
still satisfying the assumptions of Mie Theory (jBazell fe Dwek!ll990f ). 



Finally, the flux (in W cm -2 /im -1 ) from each mineral is calculated by integrating over 
the grain size distribution between the minimum radius a Q = 0.1 /im and the maximum grain 
radius a max = 100 /im, 



F em (\,r h ) 



A 2 



n(a)7ia 2 Q em (a, X)B(X, T d (a))da 



where N p is the number of particles at the peak of the grain size distribution if the peak 
of n(a) is normalized to unity, A is the geocentric distance of the dust grain, n(a) da is 
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the grain size distribution, Q e m(o>,\) is the emission efficiency for a grain of radius a(/xm) 
at wavelength A(yum), and T d (K) is the heliocentric distance dependent (rvdependent) dust 
temperature for a single grain size, a, calculated by radiative equilibrium. The flux from 
each mineral species is coadded and fit to the Gemini spectra by varying N p . The best-fit is 
determined by calculating the lowest xl (reduced \ 2 ) value at a confidence level of two-sigma. 



RESULTS & DISCUSSION 



4.1. Model Fits 



The best-fit model parameters for SW3-[B] are listed in Tableland the spectral decom- 
position is shown in Fig. El Similarly, the best fit model parameters for SW3-[C] are listed 
in Table El and the spectral decomposition is shown in Fig. |H 

The grain properties of both SW3 fragments are similar. The grains range from being 
solid to moderately fractally porous (D = 3.0 — 2.727). Fractally porous coma grains in SW3 
fragments are also required to explain the high fractional polarization measured in the optical 
I- and near-IR ij"-band a t high phase angles and throughout the comae of SW3-[B] and SW3- 
[C] by lJones et al.l ( 120081 ). The peak of the HGSD for fragment SW3-[B] in the anti-sunward 
direction is slightly larger (~ 0.5 /jm) than that of fragment SW3-[C] (~ 0.3 /xm). However, 
the grain population in SW3 overall is not dominated by extremely s mall particles as was 



evident in the comae C/1995 01 (Ha l e-Bopp) or 9P post -Deep Impact ( IWooden et al.lll999 
Harker. Woodward, fc Wooden! 1200,4 Marker et all 120071 ). 



In both fragments, most of the thermal emission is dominated by amorphous (glassy) car- 
bon and silicates. The silicate emission feature in both fragments is dominated by amorphous 
pyroxene grains, although there is some evidence for a population of amorphous olivine grains 
associated with the coma toward the nuclei of both fragments. Orthopyroxene is seen in the 
extended coma in both the sunward and anti-sunward direction of fragment SW3-[B], but 
only in the sunward direction, close to the nucleus (1") in fragment SW3-[C]. A modest pop- 
ulation of crystalline olivine grains is present in both fragments. The silicate cr ystalline mass 



fraction for the submicron to micron-size portion of the grain size distribution (IHarker et al. 
20021 iMoreno et al.ll2003l ). defined as f crys t = (crystalline) /(crystalline + amorphous), for 



SW3-[B] ranges from 0.3351K 

at the nucleus, declining to 0.061_ ' 061 at a cometocentric 



derived for SW3- 



cryst 

+ q^ 73 at a cometocentric distance of 



distance of 6.0" out into the coma in the anti-sunward direction. The /, 
[C] is 0.257to;o43 towards the nucleus, declining to 0.0 
5" into the coma in the anti-sunward direction. Therefore, within the 1-sigma confidence 
level, both fragments have a similar crystalline silicate fraction. For many EC comets and in 
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iv±ji,-iiuLi ±o ±i±£ 

Harker et al. 


11 \ .lrr 

2007; 


y.it. — u.uo u 

Sueita et al. 


2005) 



9P, fcryst — 0.3 was found in the post-impact cor 
Overall, lack of significant mineralogical differences between SW3 fragments may indicate 
that the original nucleus of t he comet was g enerally homogeneous in composition as opposed 
to that of lP/Halley or 9P jWoodeiJliool l. 



The spectra of fragment SW3-[C] in the sunward direction are dominated by emission 
from amorphous pyroxene and carbon, with little to no emission from crystals (Table ED- The 
peak of the grain size distribution is also slightly larger (0.5 -1.3 jttm) compared with the 
derived values in the anti-sunward direction (0.3 - 0.4 /mi), in addition to being more solid 
than the latter. We attribute these differences in the grain properties to the dynamics of 
the grains as they are released from the nucleus of fragment SW3-[C]. Smaller and/or more 
porous grains which can exhibit crystalline emission are more easily affected by the radiation 



pressure or are better c oupled to and entrained in the sublimating gas (e.g., (IHarker et al. 



20071 ; lOrosei et al.lll995l ) ) . Likely these grains are more easily swept out into the anti-sunward 
direction. However, a full dynamical analysis of the trajectory of the grains is outside the 
scope of this work. 

The relative abundance of amorphous carbon in the anti-sunward coma of SW3-[B] 
ranges between ~ 45 — 23 wt%, while for SW3-[C] the value ranges between ~ 60 — 42 wt%. 
The wt% of amorphous carbon in both fragments varies a little with distance from the 
nucleus (~ 20%) and is com parable to me an values deduc ed for other c omets , espe cially 
ECs (10 to 50 wt%) comae flWoodenl lioosh . For example barker et all tooi . 120041 ) find 
amorphous carbon abun dance (relative t o oth er refractory materials) of 21 wt% for comet 
C/1995 01 ( Hale-Bopp), IWooden et all J20041 ) cite 15 wt% for C/2001 Q4 (NEAT), while 
Harker et al.l (120071 ) deduce a 28 wt % relative abundance for the pre-impact coma of 9P. Each 
of the later values (including those for SW3) were derived using the same discrete composition 
thermal model approach for SED decomposi tion, enabling a consistent comparison of relative 
abundances that are not model dependent (IWooden!l2008l ). 



The average silicate-to-carbon ratio in the extended anti-sunward coma for fragment 
SW3-[B] ranges between 1.156 — 3.43, while for fragment SW3-[C] the value ranges between 
0.671 — 1.882. Overall, these ratios are comparable to yalues derived for other EC comets 
jKellev et alJbood : barker et al1l2007[ iLisse et al-lliooS iReach et al.lboioh . Fragment SW3- 
[B] went through an outburst event just prior to our observations and this event could be 
responsible for the higher silicate-to-carbon ratio in fragment SW3-[B] compared to that of 
fragment SW3-[C], as well as the slight differences in the peak of the grain size distribution 
between the two fragments. The date of the outbu rst from fragment S W3-[B] has been 
estimated to have been between 2006 April 16 to 26 (IVincent et al.ll2010l ). Likely, material 
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from the outburst is responsible for the region of enhanced flux in fragment SW3-[B] seen 3" 
from the nucleus region in our images. The activity from fragment S W3- [C] is primarily from 
two active jets on the nucleus (I Vincent et al.ll2010f ) that could be preferentially ejecting small 
grains from the nucleus subsurface. On the other hand, fragment SW3-[B] was exhibiting 
fragmentation that released all subsurface grains in a single event. Therefore, in bulk, the 
grains released during a fragmentation event are slightly larger and more silicate rich than 
those entrained in the active jet outflows from the nucleus surface. 

The radial difference in the silicate-to-carbon ratio in each fragment could arise from 
fragmentati o n of c omposite silicate rich grain aggregates entrained in the out flowing gas. 
Jones et al.l (120081 ) observed fragments SW3-[B] and SW3-[C] in imaging polarimetry mode 
in the near-IR, finding that the surface brightnesses of both fragments are consistent with 
significant grain fragmentation. Their model requires a factor of 10 change in the mean grain 
size within 2-200 km from the nucleus (corresponding to 2" in our Michelle observations). In 
addition, the scattered light polarization profile slightly increased over the same distances. 
The latter observations suggest that the grains were fragmenting, although the polarization 
could have also increased from the release of volatiles such as water ice andorganics, which 
may change the light scattering properties of the grains (j Jones et al.l 120081 ). Despite the 
evidence for fragmentation in the scattered light, the grain size distribution derived from our 
thermal models do not appreciably vary, nor does a p . Our best-fit grain size distributions 
have N = 3.4. This shallow power-law slope suggests that the scattered light could be 
dominated by the largest grains which are too cool to be seen in the 10 fim region. Because 
smaller grains have higher radiative equilibri um temperatures, the mid-IR emission at 10 /im 
may be dominated by ~ 0.1 — 1 //m grains (IReach et al.l 120071 ) . Therefore the two datasets 
may not be so incongruent if the largest grains, perhaps 10s or 100s of fim in size, are 
fragmenting, but do not significantly affect the 8-13 SED because of their cool radiative 
equilibrium temperatures. 

Alternatively, the increase in silicate-to-carbon ratio with distance from the nucleus 
could be due to grain sorting in the coma. In this scenario, the more transparent silicate 
species would dynamically separate from the dark amorphous carbon grains. Moderately 
porous (D = 2.857) amorphous carbon grains ha ve a 2-3 time s larger response to radiation 
pressure than do the same-sized silicate grains (jKelleyi 120061 ). Grain sorting by size and 
mineral composition was suggested from the time-of-fl ight analysis of the 9P post-impact 
coma contaminated by material from the ejecta plume (IHarker et al.l 120071 ). 
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4.2. Comparison to Other Comets 



It is instructive to compare the spectral signature of coma dust in the environs of the 
SW3 fragments to other comets observed at mid-IR wavelengths to place the conclusion 
based on thermal model analysis and mineral decomposition analysis in context. Fig. [5] is 
a composite plot of the spectra of fragments SW3-[B], SW3-[C], comets C/1995 01 (Hale- 
Bopp), 9P post-Deep Impact, and C/2001 Q4 (NEAT) observed with ground-based mid-IR 
spectrographs, C/2004 Bl (LINEAR) and 17P/ Holmes observed with Spitzer IRS spectro- 
graph; and lP/Halley observed by the Kupier Airborne Observatory. Comets C/1995 01 
(Hale-Bopp) and 17P/Holmes (17P) clearly have the strongest silicate features (defined as 
the excess over the blackbody continuum). However, we note that 17P was observed by 
Spitzer a little more than a couple of weeks after an extremely violent outburst event. This 
material was still within the Spitzer IRS spectral beam, and is therefore the reason for 
the relatively large silicate feature. The grains in the coma of C/1995 01 (Hale-Bopp) 
were fr act ally very porous (D = 2.5) with the HGSD peaking at a p = 0.2 /im ( Harker et al. 



20021 ). The coma dust in C/1995 01 (Hale-Bopp) is also mineralo gically diverse and includes 
amorphous and crystalline forms of both olivine and pyroxene (IHarker et al.ll2002l ). Simi- 
larly, the dust ejected from 9P post-Deep Impact was also mineralogically diverse and small 
(a p = 0.2 jum) in size, although the grains were moderately fractally porous (D = 2.727) 
Marker. Woodward, fc Wooden! 120051 ) . Furthermore, comet C/2001 Q4 (NEAT) also dis- 
played a mineralogically diverse grain population, however, the grains were solid (D = 3.0 ) 
and and slightly larger than those in Hale-Bopp and 9P (a p = 0.3 /im) (IWooden et al.ll20041 ). 



Comet 17P/Holmes also was observed to have distinct Mg-rich crystalline olivine fea- 
tures at 11.2 /im and 11.9 /im in apertures centered on the nucleus and in the diffuse coma, 



as well as toward the isolated dust cloud arising from the violent outburst event (IReach et al. 



2010l ; IWatanabe et al.ll2010l ). In 17P, dust producing these 10 /im features may arise from 



fluffy aggregate complexes containing small Mg-rich crystalline silicates (e.g.. iKimura et al. 



20081 ). or small silicates (a<0.1 /im) admixed larger amorphous dust. Sharp crystalline 
silicate features from small grains superposed on a thermal continuum arising from larger 
amorphous carbonaceous and silicate species in the 10 /im SED were evident in remote sens- 



ing o bservations of comet lP/Halley by the Kupier Airborne O bservatory (IBregman et al. 



19871 ). while the in situ flyby measurement through the coma ([McDonnell et al.l 119871 ) in- 
dicated that the grain sized distribution was dominated by large ~ 1 /im sized particles. 
Comet C/2004 Bl (LINEAR) has a relatively weak and featureless silicate feature indicative 
of emission from a large (> 1/rni) grain population. 

While there are similarities in dust grain characteristics between C/1995 01 (Hale- 
Bopp), 9P, and 17P possibly a signature of a common formation zone in the solar nebula of 
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the parent body aggregates for these comets, the composition of SW3 seems distinctly differ- 
ent. As SW3 fragmented into may pieces that subsequently sublimated leading to outbursts 
of gas entrained with dust particles from within the bulk interior of the original nucleus 
(natural excavation at depth akin to the subsurface penetration of Deep Impact projectile), 
surface weather affects cannot be invoked to explain dissimilarities in grain properties. ECs 
with frequent perihelion passages and low semi-major orbital axis are thought to have coma 
grain populations dominated by larger, com pact dust particles arisin g; from highly processed 
surface layers constituting the nucleus crust (IPrialnik fc Meklerlll99ll ) as opposed to the more 
pristine grains released in active Oort cloud comets. Larger grains are necessary to explain 
the aperture- dependent polarization properties and l ow total polarization ( < few percent) 
of some ECs (I Jones et al.ll2008l ; iKoloklova et al.ll2007l ). However, SW3 has an unusually high 
polarization for an EC in both the optical and near-IR at a range of phase angles, which lead 
Jones et all (120081 ) to suggest that observational feature is best explained by the release of 
unprocessed, porous aggregates with small effective radii from the interior of disintegrating 
comet nuclei. 

Although fragment SW3-[B] underwent an outburst event revealing pristine grains from 
the nucleus interior, the grain properties and mineralogy of its dust does not resemble that 
of either 9P or of Hale-Bopp. This dissimilarity could be attributed to 1) the smaller, more 
porous dust grains traveling out of the coma of fragment SW3-[B] prior to our observations; or 
2) the dust properties of SW3 fundamentally differ from that of 9P (and Hale-Bopp) implying 
a different formation process or evolutionary history. However, the former conjecture is 
ruled out by polar imetric observations and dust fragmentation studies of the SW3-[B] coma 
(IJones et al.ll2008l ) obtained at epochs contemporaneous with our mid-IR spectrophotometry. 



5. CONCLUSIONS 

Mid-IR images and spectroscopic observations of the ecliptic (Jupiter family) comet 
73P/Schwassman-Wachmann 3, fragments SW3-[B] and SW3-[C], were analyzed using ther- 
mal emission models. Based on our modeling, the comae dust, and by inference the bulk 
internal properties, of the individual fragments are similar with the observed mid-IR flux be- 
ing dominated by emission from amorphous pyroxene and amorphous carbon, with evidence 
for emission from amorphous olivine towards the nucleus in both fragments. Emission from 
both fragments shows evidence for the existence of some crystalline olivine in the extended 
coma in the anti-sunward direction. The grain properties in both fragments are also similar: 
solid to moderately porous grains (D = 3.0 — 2.727) with a large grain slope of N = 3.4 
and the peak of the HGSD ranging between 0.3 and ~ 1.0 /im. Finally, the two fragments 
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exhibit a similar, but slightly different silicate-to-carbon ratio (1.341 for SW3-[B] and 0.671 
for SW3-[C]); a ratio equivalent to those found for other ECs. The slight differences in grain 
size distribution and silicate-to-carbon ratio between the two fragments likely arises because 
SW3-[B] was actively fragmenting throughout its passage while the activity in SW3-[C] was 
primarily driven by jets. The breakup event populated the coma of fragment SW3-[B] with 
subsurface material. Conversely, for fragment SW3-[C], the dust population of its coma was 
driven by two active jets. 

There are other differences between the two fragments, including the detection in frag- 
ment SW3-[B] of orthopyroxene in both the sunward and anti-sunward directions. In frag- 
ment SW3-[C], orthopyroxene is primarily detected in only the sunward direction. The 
crystalline silicate fraction is higher in SW3-[B] than in SW3-[C] based on constraints de- 
rived for the crystalline silicate fraction from thermal modelling of the SEDs. Although f cryst 
for both fragments is similar to the range of fcryst inferred for other ECs. 

Comparison of coma dust properties of 73P to those of derived for ecliptic comet 
9P/Tempel 1 just after the Deep Impact event, and to the archetypal, "pristine" nearly- 
isotropic (Oort cloud) comet C/1995 01 (Hale-Bopp) suggests that materials released from 
fragment SW3-[B] which may expose grain material buried at depth below the surface of 
the original larger nucleus is devoid of very small (a < 0.5 fim) dust grains. The apparent 
absence, or diminished population of such small grain species may either be a result of dy- 
namical sorting whereby this population has passed out of the coma before our observations, 
or more likely, the initial parent body makeup of 73P is different than that of 9P/Tempel 1 
or even C/1995 01 (Hale-Bopp). 

Data discussed in this manuscript are based on observations obtained at the Gemini Ob- 
servatory, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the 
National Science Foundation (United States), the Science and Technology Facilities Coun- 
cil (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), Ministerio da Ciencia e Tecnologia (Brazil) and 
Ministerio de Ciencia, Tecnologia e Innovation Productiva (Argentina). DEH and CEW ac- 
knowledge support for this work from the National Science Foundation grant AST-0706980. 
DEH, DHW and CEW also acknowledge partial support for this work from NASA Planetary 
Astronomy Grant RTOP 344-32-21-04. MSK acknowleges support from NASA Planetary 
Astronomy Grant NNX09AF10G. The authors also would like to thank the Gemini Observa- 
tory staff for their support in conducting these observations as well as an anonymous referee 
whose comments improved the manuscript. 
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I 

to 

a Start time of integration. 10 



b Integration time on-source. 

c Low N = Low-Res 10 /im spectrum; Low Q = Low-Res 20 /im spectrum. 
d Slit rotated to position angle = 25° 
e Slit rotated to position angle = 45° 

f Flux in a 7.44" diameter aperture centered on the peak isophote of the surface brightness distribution. 



Table 2. BEST-FIT THERMAL EMISSION MODEL PARAMETERS FOR FRAGMENT B 



iV p (xl0 16 ) 

Extraction 3, Offset b N M a p c D Amorphous Amorphous Amorphous Crystalline Ortho- xt 
Aperture (arcsec) Pyroxene Olivine Carbon Olivine Pyroxene 



A +1 3.4 27.20 0.9 3.000 0.6431^?? O.OOOlg;^ O.OOOt&ooo 0.000±g;g^ - 854 ia233 2 - 71 

B 3.4 13.60 0.5 2.727 3.310±g;g|g 0.782+°;^ 9-703±°;jj^ 0.686±g:Se 2648 -L3?i °- 75 

C —I 3.4 13.60 0.5 3.000 2.671^;^ 0.000lg;ggg 5.296+^^ 0.268lg;^jj 1.701+g;^ 1.14 

D -2 3.4 10.20 0.4 2.857 5.933^;^ 0.000t{j;Sjj(j 9.800^369 0.980^;^ 1.540^;^ 1.04 

E -3 3.4 13.60 0.5 2.857 3.62]+°^ O.OOOlg;^ 4.913^223 0.170^;^ 1.481^;^ 1.55 

F -4 3.4 13.60 0.5 3.000 2.827^;^ O.OOOljj;^ 4.8371J53" 0.393±g;l?? 2.159^;^ 1.00 

G -5 3.4 13.60 0.5 3.000 3.949lg:SSa °- 000 ^o:ooo S.lOllJj;^ 0.4091°;^ 1.1451°;^ ,1.38 

H -6 3.4 13.60 0.5 3.000 4.494l£2?? 0.000±g;ggg 1.877±g;igi O.OOOlg;^ 0.384lg;^ ^.46 



a Lable of spectral extraction aperture (see Fig. [2]) 

b Distance that the center of the 0."6 x 1.0" extraction box is offset from the peak brightness. A positive value indicates an 
offset in the N-E direction, and a negative value indicates an offset in the S-W direction (see Fig. [2]). Both offsets are along the 
slit orientation axis of PA = —25 (see text). 

c Derived parameter 

d Bold values indicate values constrained to a confidence level of 1-sigma. 



Table 3. BEST-FIT THERMAL EMISSION MODEL PARAMETERS FOR FRAGMENT C 
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a Same as for Table [2] 
b Same as for Table [2] 
c Derived parameter 

d Bold values indicate values constrained to a confidence level of 1-sigma. 



Table 4. MASS OF SUBMICRON GRAINS FOR FRAGMENT B 
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a ) Defined as /c* c t ates = (crystalline)/(crystalline + amorphous). 



Table 5. MASS OF SUBMICRON GRAINS FOR FRAGMENT C 
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a ) Defined as ff icates 



' cryst 



(crystalline) /(crystalline + amorphous). 
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Fig. 1. — Images of comet 73P/Schwassmann-Wachmann 3 fragment SW3-[B] (top row) and 
fragment SW3-[C] (bottom row) at 11.6 fim (left column) and 18.1 /im (right column). The 
peak isophotes in the two panels for SW3-[B] represents 17.4 mJy, and for SW3-[C] represents 
76.6 mJy. The ratio between isophotes in all panels is a/2. The orientation of the frame on 
the sky, the direction of the Sun, and the direction of the velocity vector are indicated by 
arrows. The spatial scale of the panels for SW3-[B] is 1828 x 1828 km and for SW3-[C] is 
1697 x 1697 km. 
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Fig. 2. — Spectral extractions of comet 73P/Schwassmann-Wachmann 3 fragments SW3-[B] 
(top) and SW3-[C] (bottom). The left hand panels are 11.6 /zm images of the fragments which 
show the location of the slit and the size and location of the spectral extraction apertures 
(labeled with letters). The right hand panels show the 10 /zm spectra for the fragments. For 
SW3-[B], the spectral extractions for locations C and D have been multiplied by a factor of 
1.5 and 1.2, respectively, for display purposes. For SW3-[C], the extractions have been split 
to show extractions North-East of the nucleus (which is at extraction point M, and is shown 
on both sides for context) (left) and South- West of the nucleus (right). Also shown in the 
right hand panels are the 20 /zm spectral extractions for fragment SW3-[B] and the 18.1 /jm 
flux points for fragment SW3-[C]. 
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Fig. 3. — Spectral decomposition of model fits to the coma thermal spectra of comet 
73P/Schwassmann-Wachmann 3 fragment SW3-[B]. The location of the extraction beam 
is indicated in each panel. Overlayed on the spectra (black circles) is the total SED fit (red 
line). The minerals used to create the SED are: amorphous pyroxene (blue line), amorphous 
olivine (cyan line), amorphous carbon (orange line), crystalline olivine (green line), and or- 
thopyroxene (violet line). A solid line for a particular mineral indicates that the mineral is 
constrained by the model, a dashed line indicates that the mineral is not constrained. The 
three vertical black dashed lines indicate the wavelength locations for the resonance features 
of orthopyroxene at 9.3 and 10.5 /im and for crystalline olivine at 11.2 /im. 
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Fig. 4. — Spectral decomposition of model fits to the coma thermal spectra of comet 
73P/Schwassmann-Wachmann 3 fragment SW3-[C]. The location of the extraction beam 
is indicated in each panel. Overlayed on the spectra (black circles) is the total SED fit (red 
line). The minerals used to create the SED are: amorphous pyroxene (blue line), amorphous 
olivine (cyan line), amorphous carbon (orange line), crystalline olivine (green line), and or- 
thopyroxene (violet line). A solid line for a particular mineral indicates that the mineral is 
constrained by the model, a dashed line indicates that the mineral is not constrained. The 
three vertical black dashed lines indicate the wavelength locations for the resonance features 
of orthopyroxene at 9.3 and 10.5 /im and for crystalline olivine at 11.2 /im. 
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Fig. 5. — Blackbody normalized comet spectra extracted from slit positions coincident with 
each nucleus. The normalization is derived from fits to continuum spectral points < 8.2 /im 
and > 12.5 /im lying outside the silicate band. The blackbody normalized spectrum of 
comet 73P/Schwassmann-Wachmann 3 fragments SW3-[B] (red; B A (T) = 289 K, r^ = 1.11, 
A = 0.15) and SW3-[C] {blue; B A (T) = 298 K, i h = 1.09, A = 0.13) 10 /im spectra are 
compared to other ecliptic (Jupiter family) and nearly isotropic (Oort cloud) comets observed 
with similar ground-based remote sensing techniques (except for C/2004 Bl (LINEAR) and 
17P/Holmes which were observed with Spitzer and lP/Halley which was observed with the 
Kupier Airborne Observatory): C/2004 Bl (Linear) (open black circles; B A (T) = 215 K, 
i h = 2.06, A = 1.60), C/2001 Q4 (NEAT) (orange; B A (T) = 310 K, r h = 0.97, A = 0.35), 
9P/Tempel-1 post-Deep Impact (green; B A (T) = 288 K, i h = 1.51, A = 0.89), lP/Halley 
(open red circles; B A (T) = 300 K, r h = 1.32, A = 0.77), C/1995 01 (Hale-Bopp) (black; 
B A (T) = 288 K, r h = 0.93, A = 1.46), and 17P/Holmes (open blue circles; B A (T) = 206 K, 
r^ = 2.51, A = 1.87) The dashed lines indicate the wavelengths of crystalline silicate emission 
features at 9.3 and 10.5 /mi (orthopyroxene), and 11.2 /im (crystalline olivine). 



